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ABSTRACT 


The  seismic  surface  wave  and  P-wave  data  generated  by 
the  Bilby  explosion  and  the  associated  cavit  /  collapse  are 
studied  comparatively  to  determine  the  radiation  patterns 
of  these  waves.  The  asymmetric  radiation  patterns  of  P 
and  Rayleigh  waves  as  well  as  the  presence  of  Love  waves 
are  explained  in  terms  of  a  composite  source.  This  consists 
of  an  isotropic  dili.tational  component  due  to  the  explosion 
and  a  double-couple  component  due  to  tectonic  effects.  The 
relative  strength  of  the  multipolar  component  is  0.47  timer; 
that  of  the  explosion.  The  source  time  functions  of  Rayleigh 
and  Love  waves  from  the  Bilby  are  determined.  For  Rayleigh 
waves  this  is  a  pulse  of  the  form  p(c)=t.  exp(-1.5t_). 

For  Love  waves  the  source  time  function  may  be  a  step  function 
with  a  slow  rise  time. 


1 .  INTRODUCTION 

The  underground  nuclear  explosion  Bilbv  was  detonated  in  the 
Nevada  Test  site  on  13  September  1963  at  17:00:00  UY.  The  shot  was 
at  a  depth  of  2314  ft.  below  the  surface  and  the  coordinates  of  the 
epicenter  were  37°  03'  38"  N  and  116D  01'  18"  W.  The  geologic  medium 
at  the  shot  was  tuff.  The  assigned  magnitude  of  the  event  was  m  *  5.8 

Bilby  generated  seismic  waves  which  were  recorded  well  at  most 
of  the  North  American  stations  and  especially  at  the  LRSM  stations. 
Furthermore,  it  was  followed  by  the  cavity  collapse  (assigned  magni¬ 
tude  -  4.5)  at  17:31:  (20.5)  U  T.  The  presence  of  this  collapse, 
the  relatively  large  size  of  the  explosion,  radiation  of  both  Rayleigh 
and  Love  waves,  and  the  proximity  of  the  site  to  those  of  other  well- 
recorded  explosions  such  as  Haymaker  and  Sedan  have  made  Bilby  ideally 
suited  for  the  study  to  determine  the  mechanism  of  generation  of 
seismic  waves  and  the  radiation  pattern. 

In  this  study  we  follow  a  procedure  of  analysis  similar  to  that 
used  in  the  earlier  studies  (Tokstiz,  et  al . .  1964,  and  1965).  To 
minimize  the  effects  of  the  propagation  we  determine  radiation 
patterns  of  the  P-waves  and  Rayleigh  waves  of  Pilby  relative  to  those 
of  the  collapse.  Then  we  use  the  ratios  of  the  Love  and  Rayleigh 
wave  amplitudes  to  determine  a  source  mechanism  which  is  consistent 
with  the  Rayleigh  wave  radiation  pattern.  The  source  time  function 
is  determined  from  the  amplitude  spectra  corrected  for  propagation 
and  instrument  effects.  Finally,  we  compare  the  Bilby  results  with 
those  of  other  explosions. 

H  •  RADIATION  PATTERNS  OF  SURFACE  AND  BODY  WAVES 

The  nuclear  explosion  Bilby  generated  Love  waves  which  were 
recorded  well  at  a  number  of  the  LRSM  stations  shown  in  Figure  1. 

Since  the  long-period  horizontal  instruments  are  generally  oriented 
in  a  radial  and  transverse  direction  relative  to  the  explosion  site, 
the  Love  waves  can  be  identified  and  separated  without  much  diffi¬ 
culty.  Furthermore,  for  the  continental  paths,  the  Love  wave  group 
velocities  are  generally  higher  than  those  of  Rayleigh  waves  in  the 

period  range  of  10  -  40  sec.  This  also  facilitates  the  study  of  the 
Love  and  Rayleigh  waves. 


The  generation  of  these  transverse  waves  by  Bilby  suggests  a 
source  mechanism  which  was  more  complicated  than  a  theoretical  ex¬ 
plosive  source.  The  possible  sources  of  the  Love  waves  generated 
by  explosions  were  discussed  in  earlier  papers  (Tokaflz,  1966) . 

There  is  strong  evidence  that  the  Love  waves  are  the  results  of 
radiation  of  some  tectonic  strain  energy  either  because  of  the  ex¬ 
plosion-formed  cavity  or  the  explosicn-induced  directional  cracking 
or  faulting.  (See  Press  and  Archambeau,  1962;  Archambeau,  1964; 
ToksiJz,  et  al.,  1965;  and  Toksdz,  1966,  for  details).  Then,  the 
source  must  be  treated  as  a  combination  of  a  radially  symmetric 
explosion  and  a  multipolar  source  of  tectonic  strain  release. 

A.  Absolute  Amplitudes 

The  recorded  amplitudes  of  body  and  surface  waves  can  be 
corrected  for  the  seismograph  response  and  the  ground  displacements 
can  be  determined.  Under  ideal  conditions  these  can  be  interpreted 
in  terms  of  the  source  properties  after  corrections  for  the  geometric 
spreading  and  attenuation.  In  practice,  however,  the  geologic  con¬ 
ditions  are  complicated  and  the  idealized  situation  cannot  be  realized. 
First,  the  crustal  layering  between  the  source  and  the  station  does 
not  remain  perfectly  uniform  and  horizontal.  Second,  the  crustal 
structures  under  the  stations  affect  the  amplitudes,  and  these 
structures  could  vary  significantly  from  one  location  to  another. 

With  the  full  realization  of  these  limitations,  we  made  an 
attempt  to  test  the  possibility  of  determining  the  radiation 
pattern  from  the  absolute  amplitudes.  Figure  2  shows  the  polar 
plot  of  the  ground  motion  amplitudes  of  the  Z  component  of  Rayleigh 
waves,  if  the  radiation  pattern  and  the  structures  were  uniform, 
one  would  expect  circular  amplitude  contours  falling  off  as  ( 2rrr )** 
away  from  the  center.  The  figure  indicates  amplitude  discrepancies 
along  a  single  azimuth.  These  cannot  be  explained  by  the  source 
properties  and  must  be  the  results  of  the  structures  and  ground 
coupling  effects  at  the  stations. 

Another  example  of  these  structure  effects  are  illustrated 
in  Figure  3  where  magnitudes  (m)  are  plotted  against  the  azimuth 
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in  a  polar  diagram.  The  data  are  from  the  Bilby  shot  report, 

(Teledyne,  1963) .  Again  the  scatter  is  such  that  no  radiation 
properties  can  be  determined.  In  the  case  of  body  waves,  the  mantle 
as  well  as  the  crustal  structure  plays  a  significant  role  in  con¬ 
trolling  the  amplitudes.  At  epicentral  distances  corresponding  to 
the  depths  of  velocity  discontinuities  in  the  mantle,  these  effects 
are  especially  important  since  the  amplitude  is  a  strong  function  of 
the  velocity  gradient  (dv/dr)  at  depth  (Asbel,  et  al . .  1966;  Chinnery 
and  Toksdz,  1966;  Anderson,  1966).  Our  knowledge  of  the  exact  mantle 
structure  along  each  path  is  very  poor,  and,  at  the  present,  we 
cannot  hope  to  compute  these  effects  to  make  the  necessary  corrections. 
B.  Relative  Amplitudes  of  Rayleigh  WaveB. 

Since  the  effectiveness  of  using  the  absolute  amplitudes  of 
seismic  body  and  surface  waves  are  limited  by  the  crust  and  mantle 
structures,  a  method  of  normalization  to  minimize  '-.hese  effects 
is  considered.  Unfortunately,  the  source  and  the  propagation  factors 
cannot  be  isolated  from  each  other.  Thus,  we  must  consider  methods 
in  which  one  of  the  variables  is  held  constant  while  the  effect  of 
the  other  is  investigated.  In  this  work  we  uBed  two  methods  for 
the  study  of  relative  amplitudes  tc  determine  the  radiation  pattern 
from  Bilby.  In  the  first,  we  normalized  the  amplitudes  of  ..he 
Rayleigh  and  P-waves  generated  by  the  explosion  to  those  of  the 
collapse  event  that  followed.  In  the  Becond  we  UBed  the  ratio  of 
the  Love  and  Rayleigh  wave  amplitudes  to  obtain  a  source  function. 

The  collapse  of  the  cavity  (formed  by  the  explosion)  seems  to 
provide  an  excellent  reference  for  normalization.  When  the  ampli¬ 
tude  ratios  of  the  explosion-and  collapse-generated  waves  are 
taken,  the  propagation  and  instrument  effects  completely  cancel 
out,  and  the  ratio  directly  reflects  tie  source  effects.  In  pre¬ 
vious  studies,  it  was  found  that  ti.a  radiation  pattern  from  the  post¬ 
explosion  collapse  was  in  general  more  symmetric.  The  uniform 
amplitudes  as  well  as  the  absence  of  outstanding  Love  waves  were 
the  results  of  this  radial  symmetry  (ToksOz,  et  al..  1965;  Toksdz, 

1966)  . 
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If  we  assume  that  this  was  the  case  for  Bilby  (i.e.,  that  the 
Bilby  collapse  had  a  radially  symmetric  radiation  pattern),  then 
the  explosion/collapse  ratio  should  show  whether  the  explosion  was 
symmetric.  The  ratios  of  peak  amplitudes  of  the  Rayleigh  waves 
are  shown  in  Figure  4  as  a  function  of  azimuth.  The  peak  ampli¬ 
tudes  were  directly  read  from  the  long-period  records  of  the  LRSM 
stations.  In  a  few  instances  the  spectral  ratios  were  computed, 
and  these  in  the  average  gave  the  same  results  as  the  peak  ampli¬ 
tudes  . 

The  azimuthal  coverage  in  Figure  4  is  far  from  being  complete, 
but  in  the  north  and  northeast  directions  where  there  is  a  sufficient 
number  of  observations,  the  deviation  from  a  uniform  ratio  is  clear. 
If  we  assume  that  the  radial  non-uniformity  of  the  radiation  pattern 
from  the  explosion  was  due  to  some  form  of  tectonic  complication 
(such  as  the  relaxation  of  the  medium  due  to  the  cavity  or  induced 
rupture-  Press  and  Archambeau,  1962;  Archambeau,  1964;  ToksOz,  et  al, 
1965)  we  can  include  this  effect  by  superimposing  a  multipolar  term 
to  the  explosive  source.  Both  seismic  model  experiments  and 
theoretical  studies  indicate  that  at  large  distances  from  the 
source,  a  double-couple  type  source  function  is  a  good  representa¬ 
tion  for  tectonic  strain  release  (Honda,  1962) 

Then  the  displacements  observed  at  a  distant  station  can  be  written 
as  the  vectorial  sums  of  those  due  to  an  explosive  source  and  to  a 
double-couple . 

Using  the  notation  of  Toksflz,  et  al . ,  (1965),  we  can  write 
the  far-field  expressions  for  the  Rayleigh  wave  ground  displace¬ 
ments  from  a  r.ear-surface  explosive  source. 

C  U 

*  (ui)  -  - --i—  kY-J? — )  A_(w)  T(m)  exp(-  Y_r)  exp  [i  (wt-k  r-<p  +3n/4)  ] 

*  (2nr)S  *0  '  R  R  R  t 


f  °  0 


-J  2ARt»>)  T  (m)  exp(  -VRr)  expZi  (u>t-kRr-rt3n/4)  ] 


0 


W  (w)  ,  U  (w) ,  V  (w)  are  the  vertical,  radial  and  tangential 
e  e  e 

components  of  the  displacement,  k  is  the  wave  number,  r  is  the 

R 

radial  distance,  yd  is  the  Rayleigh  wave  attenuation  coefficient. 

R 

A  is  the  medium  response  for  Rayleigh  waves  due  to  a  vertical 
R 

force,  u  and  w  are  the  components  of  particle  velocity  at  the 
o  o 

surface.  T('i>)  and  qp  (•.«)  are  the  amplitude  and  the  phase  spectra 
of  the  source  time  function.  The  displacements  due  to  an  ortho¬ 
gonal,  horizontal  double-couple  source  are  (Ben-Menahem  and  Hark- 
rider  1964;  Toksdz,  et  al ■ ,  1965): 

* 

C  r  u 

W  (id)  -  - r —  ^ - J  Ar(ib)  sin20  exp(-  yRr)  eypfi (u>t-kRr-cp^  +  3n/4)  , 

C  (2irr)’  w 


U.  (id)  “  — - — r-  A  <uu)  T*  (id)  sir>2©  exp  (-  Y_r)  exp[i(i»t-k  r-cp’  -3n/4)]  (2) 

dc  (2nr)’  Vw0  R  R  Rt 


vdcC«) 


(2nr) 


’j 

\  (•)  T'(ei)  coi.20  exp(-  yt  r )  expf  i  (u>t-kT  r-cp'  -3tt/4)] 

L  L  ' 


The  subscript  R  and  L  refer  to  Rayleigh  and  Love  waves,  respect¬ 
ively.  The  angle  9  is  measured  counter-clockwise  from  the  principal 
plane  (i.  e.  fault  plane)  of  the  double-couple. 

The  far-field  displacements  of  Rayleigh  and  Love  waves  from  a 
composite  source  consisting  of  an  explosion  and  a  double-couple  can 
be  written  from  (1)  and  (2) . 


•j  -  k  <«)  +  w, Jib) 

Rz  e  dc 


f  _ ) 

»  ( “>  >  +  F  T  —  sin20  expfi  (f4>J  Ij 


uL  -  V<jcM 


(3) 


F  is  the  relative  strength  cf  the  double-couple,  =  rp^-cp^  is  the 
phase  difference  of  two  time  functions.  The  term  with  the  factor 
sin  29  gives  the  azimuthal  dependence  of  the  Rayleigh  wave  radiation. 
If  we  assume  that  the  source  time  functions  are  approximately 


release  (i.e.  T=T‘) 


the  same  for  an  explosion  and  the  tectonic  strain 
eg.  (3)  becomes 

U. 


Rz 


U. 


W  (1  +  F  sin2Q) 
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(4) 
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dc 


The  motion  from,  the  collapse  of  the  cavity  can  be  represented 

b*  (»Rl)coll.P-,-C3WeeX?(i,Pc>  (5) 

where  C  is  the  relative  amplitude  and  *c  is  the  phase  of  time 
function  relative  to  the  explosion.  From  previous  studies  of 
explosion  and  collapse  pairs  it  was  found  that,  for  long  periods, 
rp  „  rr  (Brune  and  Pomeroy,  1963;  Smith,  1963;  Toksflz,  et_al,  1964). 
This  result  means  essentially  that  the  collapse  can  be  considered 
an  implosion  when  we  are  dealing  with  long  period  data.  With  the 
above  formulations  and  assumptions,  the  theoretical  Rayleigh  wave 
displacement  ratios  for  a  given  explosion-collapse  pair  can  be 
written  as 


'  u„  \ 
V  R ) 


(uO 


explos lun_  _  c>ri  +  F  sin2e1 


(6) 


collapse 


In  Figure  4  the  theoretical  curve  is  computed  using  (6);  choosing 
C’,  F,  and  the  orientation  of  the  double-couple  to  fit  the  data 
best,  we  find  C'  -  12,  F  -  0.47,  and  the  reference  direction  for  the 
double-couple  principal  plane  (i.e.,  fault  plane  )  0  340  .  With 

these  data  the  fit  can  be  considered  to  be  good.  The  above  figures 
mean  that  the  explosion-generated  surface  waves  were  12  times  larger 
than  those  of  the  collapse,  and  that  the  relative  strength  of  the 
double-couple  force  was  0.47  times  that  of  the  explosion. 

The  p-wave  radiation  pattern  should  be  similar  to  that  of  the 
Rayleigh  waves.  Unfortunately  the  P-wave  data  from  the  collapse  are 
scarce,  and  the  amplitudes  not  very  reliable.  The  available  data 
(amplitude  ratios  of  explosion  to  collapse)  are  shown  in  Figure  5. 

It  is  obvious  that  the  radiation  pattern  of  Rayleigh  waves  shown  in 


Figure  4  does  not  fit  the  data:  the  amplitudes  are  too  large,  although 
the  shapes  of  the  two  curves  are  similar.  A  new  curve  based  on 
equation  (6)  is  computed  keeping  all  the  parameters  except  C1 
the  same.  A  value  of  C1  =4.8  seems  to  fit  the  P-wave  radiation 
pattern.  This  implies  that  the  amplitudes  of  Rayleigh  wave- 
relative  to  the  P-waves  was  larger  for  t  le  collapse  than  f  >r  the 
explosion.  This  can  be  explained  using  general  observations  about 
source  functions.  The  peak  of  the  source  spectrum  shifts  to  lower 
frequencies  with  increasing  size  of  the  event,  whether  it  is  an 
explosion  or  an  earthquake,  (Toksfiz,  et  1964).  it  has  also 

been  observed  that  in  short-period  recordings  the  collapse  event 
seems  to  have  a  source  spectrum  shifted  to  lower  frequencies  com¬ 
pared  to  those  of  the  explosions  of  equivalent  size  (Smith,  1963). 
However,  the  size  effect  seems  to  dominate  in  this  case,  and  the 
collapse  appears  to  be  richer  in  high  frequencies.  When  we  take 
collapse/explosion  source  spectral  ratios  over  a  wide  frequency 
range,  we  would  expect  a  frequency  dependence  unless  the  sizes 
are  comparable.  This  frequency  dependence  does  not  affect  our 
assumptions  and  calculations  when  we  are  limited  to  a  narrow  fre¬ 
quency  band . 

C .  Amplitude  Ratios  of  the  Love  and  Rayleigh  Waves 

The  Love  waves  generated  by  the  explosions  can  also  be  used  to 
determine  the  nature  of  the  source  mechanism.  According  to  our 
formulation,  the  Love  waves  will  be  due  to  the  contribution  of  the 
tectonic  component  of  the  source.  The  explosion  itse.U  will  not 
generate  Love  waves  under  ideal  conditions. 

In  determining  the  radiation  pattern,  we  will  again  use  a 
normalization  scheme  to  minimize  the  effects  of  propagation  and 
recording  uncertainities .  Since  we  do  not  have  a  pure  Love  wave 
source  that  can  be  used  as  a  reference,  we  will  normalize  to 
Rayleigh  waves.  This  scheme  was  successfully  applied  to  explosions 
and  earthquakes  by  Toksdz,  et  al .  (1965). 

The  ratio  of  the  Love  wave  amplitude  to  the  Z  component  of  the 
Rayleigh  waves  generated  by  the  explosion  can  be  written  using 
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equations  (1),  (2),  (3),  and  (4). 


f  k.' 


K.  cos20 

L _ 


'Rr 


(1  +  F  sin2G)  kjj  Ar  (Uq/Wq) 


1  expT-r ( v 


V5 


(7) 


In  writing  (7)  it  was  assumed  that  the  source  time  functions  T(w) 
and  T'(w)  were  the  same  for  both  the  explosion  and  the  double-couple 
component.  Furthermore,  as  in  the  previous  section,  it  is  assumed 
that  the  tectonic  contribution  can  be  represented  as  a  double-couple. 
Thus  the  Love  waves  are  those  generated  by  a  double-couple,  and  the 
Rayleigh  waves  are  the  vectorial  sum  of  those  due  to  the  explosion 
and  due  to  the  double-couple . 

In  using  (7)  we  must  compute  kL<  A^ ,  kR,  AR.  (ji*  ,w^)  * VR • 

These  quantities  are  functions  of  the  frequency  for  a  given  structure. 

Fortunately,  kL  AjVkR  AR  is  nearly  unity  for  the  frequency  range  of 

our  interest.  Thus  the  effect  of  the  structure  is  minimized  by  the 

normalization  process,  and  an  imprecise  knowledge  of  structure  does 

not  limit  the  applicability  of  the  method.  We  took  one  average 

structure  for  the  Western  United  States  given  by  Alexander  (1963) 

and  computed  the  amplitude  response  for  the  Rayleigh  (A  )  and  the 

R 

Love  (A  )  waves.  The  method  and  programs  of  Harkrider  (1964)  were 

Li 

used  in  these  computations.  The  results  are  shown  in  Figures  6  and 
7.  kL.kR.  and  (u*/wj  are  parameters  which  are  computed  by  all 
standard  dispersion  programs. 

In  computing  the  Love/Rayleigh  amplitude  ratios  from  the  long 
period  recordings  cn  the  LRSM  stations  we  used  the  peak  amplitudes. 

In  cases  where  spectra  were  computed,  they  peaked  at  about  T  «*  18 
seconds,  and  the  ratio  was  nearly  constant  in  the  period  range  of 
10  to  30  seconds.  These  sptctral  ratios  are  shown  in  Figure  8  for 
four  stations,  together  with  «-he  ratio  of  peak  amplitudes.  In  the 
average  the  agreement  is  good  enough  to  justify  the  use  of  peak 
amplitudes . 

All  the  available  U /u  data  from  Bilby  are  shown  in  Figure  9, 

Li  hZ 

as  well  as  the  theoretical  curve  based  on  equation  (7).  The  refer- 


ence  plane  for  double-couple  orientation  (plane  of  0  “  340°)  and 

the  relative  strength  F  -  C.47  are  the  same  values  that  were  determined 
from  the  explosion/collapse  Rayleigh  wave  ratios.  The  agreement  be¬ 
tween  the  observed  and  the  theoretical  can  be  considered  good. 

The  consistency  of  one  single  source  model  for  both  the  Rayleigh 
and  the  Love  wave  radiation  patterns  is  encouraging.  Although  this 
is  not  a  definite  proof,  it  gives  support  to  our  method  of  synthe¬ 
sizing  the  source  and  to  our  assumptions. 

III.  SOURCE  TIME  FUNCTION  OF  BILBY 

The  source  time  function  of  the  explosion  can  be  determined  from 
the  recordings  of  the  motion  at  distant  stations  by  correcting  for 
the  instrument  response  and  the  response  of  the  propagation  medium. 

For  surface  waves,  we  can  compute  both  amplitude  (including  atten¬ 
uation)  and  phase  response  if  the  structure  is  known,  in  this  study 
we  will  use  an  average  structure  and  use  only  the  amplitude  spectra, 
since  the  accuracy  of  phase  spectra  depends  very  strongly  on  exact 
knowledge  of  the  structure . 

The  Rayleigh  waves  recorded  at  three  stations  are  first  correct¬ 
ed  for  the  instrument  response  to  obtain  the  true  ground  displacement. 
These  stations  are  Kanab,  Utah,  Campo,  California,  and  Winnemucca, 

Nevada;  they  represent  excellent  azimuthal  coverage  at  fairly  close 
distances.  The  filtered  Rayleigh  wave  pulses  from  two  of  these  are 
shown  in  Figure  10  and  their  Fourier  amplitude  spectra  in  Figure  11. 
Ground  displacement  spectra  (Figure  12)  are  obtained  by  correcting 
the  spectra  of  Figure  11  for  the  instrument  response  at  each  station. 
They  represent  the  product  of  the  source  amplitude  spectrum  and  the 
response  of  the  layered  medium  (i.e.  propagation  path)  to  Rayleigh 

waves.  The  medium  response  includes  the  attenuation  effect  and  the 
source  depth. 

To  determine  the  source  function  we  must  correct  the  ground 
displacement  for  the  propagation  factor.  This  was  done  using  the 
impulse  response  of  the  medium  (Figure  5)  as  shown  in  equation  1. 


The  effect  of  attenuation  was  removed  using  y  ■  nf/UQ  where  U  is 

R 

the  group  velocity  and  Q  was  assumed  to  be  100,  independent  of 
frequency.  The  corrected  spectra  are  shown  in  Figure  13.  They 
represent  the  spectrum  of  the  source  pressure  function. 

The  interpretation  of  amplitude  spectra  in  terms  of  a  time 
function  cannot  be  done  unless  we  can  incorporate  either  the  phase 
information  or  some  other  constraint.  We  will  assume  that  the 
pressure  pulse  has  the  foxrm  p(t)  =  P  t  exp(-ht).  This  formulation 
was  discussed  in  an  earlier  study  (Toksdz,  et  al.,  1964).  The 
whole  problem  now  consists  of  determining  parameter  h  from  the  spectra. 
A  value  of  tj  “  L.5  seems  to  agree  well  with  the  observations  as  shown 
in  Figure  13.  The  corresponding  time  function  is  given  in  Figure  14. 

We  must  note  here  that  the  p(t)  we  determined  represents  the 
stress  wave  form  not  at  the  source  but  some  distance  from  the  point 
of  detonation.  Since  we  used  a  linear  theory  based  on  infinitesimal 
strains  to  correct  for  propagation  and  attenuation  effects,  our 
corrections  are  valid  only  to  the  boundary  of  the  region  where  these 
conditions  are  met.  This  may  be  at  a  distance  of  several  hundred 
or  a  few  thousand  meters  from  the  source  point.  We  must  clarify  one 
other  aspect  of  Figure  10:  the  pressure  pulse  was  based  on  data  in 
the  period  range  of  10  to  40  seconds,  so  we  could  not  see  an/  of  the 
fine  features  of  the  pulse  that  would  be  observable  primarily  in  the 
high  frequency  components. 

It  would  be  of  interest  to  determine  the  source  function  of  the 
Love  waves  generated  by  Bilby.  Choosing  four  stations  where  Love 
waves  are  well  separated  from  the  Rayleigh  wave  interference,  we 
followed  the  same  procedure  as  for  Rayleigh  waves  for  correcting 
the  spectra  for  the  instrument  response,  and  the  response  of  the 
layered  medium  for  an  orthogonal  double-couple  source.  The  resultant 
spectra  of  source  time  function  are  shown  in  Figure  15,  and  their 
shapes  are  consistent  for  all  four  stations. 

From  the  comparison  of  Figures  13  and  15  it  is  obvious  that 
source  time  function  of  the  explosion-generated  Rayleigh  waves  and 
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e  waves  are  quite  different.  The  Love  waves  seem  to  be  rich¬ 
er  n  low  frequency  components.  From  the  available  data  in  a  relative¬ 
ly  narrow  frequency  band  we  «>.not  determine  the  time  function,  it 
appears  to  be  between  a  step  function  and  a  ramp.  For  the  former, 
the  spectrum  would  be  a  linearly  increasing  function  of  the  period, 
and  for  the  latter  it  would  increase  as  T2.  a  step  function  with  a 
slow  rise  time  may  be  adopted.  The  source  spectra  of  the  Rayleigh 
and  Love  waves  generated  by  explosions  are  important  because  they 
contain  the  information  about  the  mechanism  of  the  Love  wave  generat¬ 
ion.  The  differences  demonstrated  above  are  significant.  They 
indicate  that  the  explosion  may  trigger  the  Love  wave  radiation,  but 
it  does  not  control  its  time  history. 

IV.  CONCLUSIONS 

Ir.  thl.  study  we  determined  the  .ounce  properties  of  the  Bilby 
explosion  from  the  radiation  pattern  of  a.ia„lc  waves.  „e  followed 
the  name  procdure.  that  were  u,ed  in  our  earlier  atndi.a  of  the 
Hardhat.  Haymaker.  Sedan  and  Shoal  exploaiona.  Here  we  will  e„„p,re 
the  Bilby  results  with  some  of  the  earlier  results. 

"°St  «•»«  that  Bilby.  like  the  Haymaker 

nd  Sho,  explosions.  ,e„.r.ted  some  Love  waves.  The  source  mechan- 
tn  all  case,  can  be  explained  in  term,  of  an  ideal  e*pl„.i„. 

^rc.  superimposed  over  a  tectonic  enure,  o,  double-couple  form. 

The  -  ..  .t  o,  and  relative  ,trength  of  the  double-couple  seem  to 
controlled  by  the  properties  of  the  medium  ,„d  the  orientation 
the  tectonic  axes.  Bilby  <i„  tuff,  ,„d  Haymaker  ,i„ 
were  located  about  5  km  apart.  The  radiation  pattern,  of  Hayl.i  „ 
wave,  are  almost  identical,  in  both  ca.es  the  principal  plan,  o,  the 

zjztj:  °tientea  in  ihe  airecti'>n  e  ■  34°°-  F°r  «•  s*°" 

xplosion,  which  was  fired  in  a  completely  different  area,  the  or- 
ientation  of  the  double-couple  was  in  verv  a  a 

of  *****  v  P  very  good  agreement  with  those 

thguakee  xn  the  area  (Toksdz,  e^l . ,  i965).  These  fact- 

suggest  that  multipolar  contributions  to  the  radiation  patterns  are 
controlled  by  the  ,e„er,l  tectonic  feature,  ov  the  region 
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The  relative  strength  of  the  tectonic  (double-couple)  contri¬ 
bution  to  the  radiation  pattern  seems  to  be  controlled  by  the 
pro^'erties  of  the  medium  in  which  the  explosion  is  detonated.  This 
can  be  justified  by  the  fact  that  for  explosions  in  salt  domes 
(Gnome  and  Salmon)  and  loose  alluvium  (Sedan)  the  source  functions 
did  not  have  multipolar  components  (i.e.,  F  *  0) .  For  Shoal, 
which  was  buried  deep  in  alluvium,  a  value  F  =  0.33  was  determined. 
For  Bilby  (in  tuff)  F  was  0.47,  and  for  Haymaker,  which  was  fired  in 
granite,  F  was  equal  to  0.9.  These  indicate  an  increase  of  F  with 
increasing  rigidity  and  shear  strati,  energy  capacity  of  the  medium. 
From  these  examples  we  may  be  able  to  conclude  that  the  multipolar 
component  of  the  seismic  energy  radiation  is  due  to  release  of  some 
of  the  strain  energy  accumulated  in  the  medium.  This  may  be  due  to 
the  relaxation  because  of  the  cavity  formation,  and/or  the  cracks 
formed  in  the  medium. 
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FIGURE  CAPTIONS 

Fig*  1*  Distribution  of  the  LRSM  stations  which  received  Bilby 

signals . 

Fig*  2.  Observed  ground  motion  amplitudes  (mu/sec)  of  the 

vertical  component  of  Bilby  Rayleigh  waves  as  a 
function  of  azimuth  and  distance  from  the  source. 

Fig*  3.  Polar  plot  of  the  Bilby  magnitudes  as  determined  at 

vatious  stations  (Radial  scale  is  magnitude  m) . 

Fig.  4.  The  radiation  pattern  of  normalized  Rayleigh  waves 

for  Bilby.  The  points  are  the  ratio  of  the  amplitudes 
of  explosion-generated  to  collapse-generated  Rayleigh 
waves.  The  curve  is  the  theoretical  radiation  pattern 
for  a  composite  source  consisting  of  an  explosion  and 
an  orthogonal  double-couple. 

Fig.  5.  The  radiation  pattern  of  P  waves  from  Bilby.  Points 

are  normalized  amplitudes  (ratio  of  explosion  to 
collapse).  Theoretical  curve  is  the  radiation  pattern 
of  a  composite  (explosive  plus  double-couple)  source. 
Outer  curve  is  same  as  that  of  Fig.  4,  included  for 
comparison . 

Fig.  6.  Rayleigh  wave  amplitude  response  of  the  layered  medium 

to  an  explosive  source  near  the  surface,  with  impul¬ 
sive  time  function. 

Fig.  7.  Love  wave  amplitude  response  of  the  medium  to  a  double¬ 

couple  source  near  the  surface. 

Fig.  8.  Love  wave  to  Rayleigh  wave  vertical  component  Fourier 

amplitude  spectral  ratio  at  four  stations.  Dashed  lines 
are  the  ratios  of  peak  amplitudes  measured  in  time  domain. 

Fig*  9*  Amplitude  ratios  of  the  explosion-generated  Love  and 

Rayleigh  waves  U^U  ) .  The  theoretical  curve  is  for 
the  composite  source 'described  in  Fig.  4. 

Fig.  10.  Rayluigh  wave  pulses  at  two  stations  plotted  from  digit¬ 
ized  data. 
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Fig.  11.  Fourier  amplitudes  spectra  of  the  pulses  shown  in 
Fig.  10. 

Fig.  12.  Ground  displacement  spectra  obtained  from  those  of 
Fig.  10  after  correction  for  instrument  response 

Fig.  13.  Amplitude  spectra  of  source  time  function  after 

correction  for  propagation  effects.  Circles  indicate 
data  ur.corrected  for  attenuation?  triangles  indicate 
the  data  corrected  taking  Qm100. 

Fig.  14.  The  Bilby  pressure  function  at  the  boundary  of  the 
linear  zone. 


Fig.  15.  Amplitude  spectra  of  source  time  function  for  Love 
waves  after  correction  for  propagation  effects. 
Circles  are  uncorrected;  triangles  are  corrected 
for  attenuation. 
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figure  4 


The  radiation  pattern  of  normalized  Rayleigh  waves 
for  Bilby.  The  point*  are  the  ratio  of  the  ampli¬ 
tudes  of  explosion-generated  to  collapse-generated 
Rayleigh  waves.  The  curve  is  the  theoretical  rad¬ 
iation  pattern  for  a  composite  source  consisting 
of  an  explosion  and  an  orthogonal  double-couple. 
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Figure  5 


Th«  radiation  pattern  of  P  wave*  from  Bilby 
Points  are  normalized  amplitudes  (ratio  of 
explosion  to  collapse) .  Theoretical  curve 
is  the  radiation  pattern  of  a  composite  (ex¬ 
plosive  plus  double-couple)  source.  Outer 
curve  is  same  as  that  of  Figure  4,  included 
for  comparison. 
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Figure  8.  Love  wave  to  Rayleigh  wave  vertical  component  (0_/URz)  Fourier  amplitude  spectral  ratio  at 
four  stations.  Dashed  lines  are  the  ratios  of  peak  amplitudes  measured  in  tima  domain. 
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Figure  12.  Ground  displacement  spectra  obtained  from  those 
of  Figure  10  after  correction  for  instrument 
response . 
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Figure  15.  Amplitude  spectra  of  source  time  function  for 

Love  waves  after  correction  for  propagation  effects 
Circles  are  uncorrected;  triangles  are  corrected 
for  attenuation. 
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